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Since  the  T970's,  Geometrical  Theory  of  Diffraction  fGTD)  solutions 
have  been  applied  to  compute  the  radiation  patterns  of  aircraft  antennas 
[1-7].  These  solutions  have  been  used  along  with  the  principle  of  reci- 
procity to  analyze  the  surface  current  and  charge  density  induced  on  air- 
craft [8-9].  These  previously  developed  solutions  were  employed  to  compute 
surface  current  and  charge  densities  for  aircraft  in  a free  space  environ- 
ment; however,  there  is  a need  for  this  information  when  the  aircraft  is 
closely  coupled  with  other  structures.  For  example,  the  surface  current 
and  charge  density  induced  by  a plane  wave  incident  on  a parked  aircraft. 

The  object  of  this  study  is  to  analyze  the  surface  current  and  charge 
density  induced  on  a finite  cylinder  situated  above  a perfectly  conducting 
qround  Diane  due  to  an  incident  plane  wave.  This  qeometrv  can  be  used 
to  investigate  the  mechanisms  involved  in  the  parked  aircraft  problem 
and  other  similar  problems. 


The  basic  approach  taken  in  this  study  is  to  employ  the  Geometrical 
Theory  of  Diffraction  fGTD).  The  GTD  requires  that  the  object  be  large 
in  terms  of  wavelength  such  that  various  scattering  centers  are  separated 
by  at  least  a quarter  wavelength.  Therefore  the  solutions  for  surface 
current  and  charge  densities  are  applicable  for  freguencies  down  to  around 
50  MHz  for  larger  aircraft  structures.  The  induced  surface  current  den- 

-A  * 

sity  is  given  by  J = 2n  x H where  n is  a unit  normal  to  the  surface  and 

-A-j 

H is  the  incident  magnetic  field ‘Intensity  on  the  cylinder  surface.  The 
induced  surface  charge  density  Is  given  by  p = ?en  x E where  e is  the 
permittivity  of  the  surroundinq  medium  and  E‘  is  the  incident  electric 
field  intensity  on  the  surface.  Usinq  reciprocity,  the  electric  surface 
current  on  a perfectly  conducting  aircraft  surface  can  be  determined  from 
the  radiation  patterns  of  two  mutually  orthoqonal  magnetic  current  moments 
which  are  located  tanqentiallv  on  the  aircraft  surface  at  the  point  where 
the  current  is  desired  [8].  The  electric  charqe  density  induced  cn  the 
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aircraft  can  be  similarly  obtained  from  the  radiation  pattern  of  an  in- 
finitesimal electric  current  moment  located  along  the  normal  to  the  air- 
craft surface  at  the  desired  observation  point. 


The  Geometrical  Theory  of  Diffraction  solution  used  in  this  study 
is  described  in  Chapter  II.  In  Chapter  III  the  two  dimensional  waveguide 
problem  is  considered.  The  GTD  solution  is  compared  with  a modal  type 
solution  to  show  where  the  GTD  ray  picture  is  valid.  Chapter  IV  describes 
the  principal  plane  scattering  from  a cylinder  in  the  presence  of  a per- 
fectly conducting  ground  plane. 

Chapter  V presents  a technigue  for  partially  correcting  caustics  which 
occur  in  the  solution  Dresented  in  Chapter  IV.  A summary  of  the  present 
study  and  discussion  is  given  in  Chapter  VI. 


II.  THEORETICAL  BACKGROUND 


A.  Introduction 

The  basic  technique  used  in  this  study  is  the  geometrical  theory  of 
diffraction  (GTD).  The  canonical  wedge  diffraction  geometry  is  illustrated 
in  Figure  1.  In  the  format  of  GTD,  the  total  electric  field  is  given  by 

-H  -*r  -M 

E = E + E + E . 

The  field  l^1  is  the  electric  field  directly  radiated  by  source,  the  field 
lfr  is  the  electric  field  reflected  from  the  wedqe,  and  the  field  is 
the  electric  field  diffracted  from  the  edge  of  the  wedge. 


OBSERVATION 


LIT  REGION 


Figure  1.  The  basic  GTD  wedge  diffraction  problem. 

The  structures  used  in  this  report  are  assumed  to  he  perfectly  con 
ducting  and  the  surrounding  medium  is  entirely  free  soace.  All  fields 
carry  an  exploit!  time  dependence  which  is  assumed  and  suppressed. 


B.  Geometrical  Optics  Fields 

The  incident  electric  field,  E1,  can  be  generated  by  an  arbitrary 
electric  or  magnetic  source  producing  plane,  cylindrical  or  spherical  waves. 
In  this  study,  electric  or  magnetic  current  moments  producing  spherical 
waves  are  used  along  with  the  reciprocity  theorem  to  model  surface  cur- 
rent and  charge  density.  The  far  field  pattern  of  the  source  is  written 


eVa,<J>)  = [flFIe,<j,l  + 4>Glfi,<t>)j 


where  F fe ,d> ) and  Gfe,d>)  are  source  pattern  functions  and  s is  the  distance 
from  the  source  to  the  observation  point. 

The  reflected  electric  field  from  a perfectly  conducting  curves  sur- 
face, as  shown  in  Figure  2,  is  expressed  as 


EJs)  = eVqJ  • E 


(pJ+sHpJ+s) 


where  E^Q  ) is  the  incident  electric  field  at  the  point  of  reflection 
0r.  The  dyadic  reflection  coefficient  17  is  given  by 

R = e'  er  - e e . 14) 

H ii  ii 


The  unit  vector  ex  is  the  vector  perpendicular  to  the  plane  of  incidence 
and  the  vectors  e,|  and  e[[  are  the  incident  and  reflected  unit  vectors  par- 
allel to  the  plane  of  incidence.  The  plane  of  incidence  is  defined  as 
the  plane  containinq  the  incident  ray  and  the  normal  Ini  to  the  surface 
at  the  point  of  reflection.  The  point  of  reflection,  0p,  is  located  using 
the  laws  of  reflection  which  state  that  the  angle  of  incidence  equals  the 
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Figure  ?.  Reflection  by  a smooth  convex  surface. 

angle  of  reflection  and  also  the  incident,  reflected,  and  surface  normals 
are  coplanar.  These  laws  can  be  expressed  mathematically  as 

- 7 • n = s • n (F) 

and 

A A A A 

I x n = s x n . ffil 

The  above  unit  vectors  are  illustrated  in  Fiqure  ?.  The  parameters  p!^ 
and  p^  are  the  princiDal  radii  of  curvature  of  the  reflected  wavefront 
at  the  reflection  Doint.  Kouvoum.iian  [in]  qives  a procedure  for  findinq 
the  radii  of  curvature  bv  diagonal izinq  the  curvature  matrix  of  Deschamps 
pi  I for  the  reflected  wavefront.  Let  the  wavefront  be  incident  on  a curved 

A A 

surface  at  the  point  Qr  as  shown  in  Figure  ?.  The  vectors  e1  and  e?  are 
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the  unit  vectors  in  the  principal  directions  of  the  surface  at  Qr  with 
r^  and  r^  being  the  two  principal  radii  of  curvature.  The  principal  radii 
of  curvature  of  the  reflected  wavefront  as  given  by  Kouyoumjian  [10]  are 
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An  asymptotic  solution  for  the  diffraction  from  a conducting  wedge 
was  first  solved  by  Sommerfeld  [12].  Originally  plane  wave  diffraction 
coefficients  as  presented  by  Keller  [13]  were  used  as  the  basis  for  the 
GTO  solution.  However,  as  shown  in  Reference  [14],  the  use  of  cylindrical 
wave  diffraction  is  necessary  for  antenna  applications.  Thus,  different 
formulations  for  wedge  diffraction  were  substituted  for  the  plane  wave 
diffraction  coefficient  which  is  the  basis  for  wedge  diffraction  theory. 
Pauli  [IF]  introduced  the  Vg  function  as  a practical  formulation  for  a 
finite  anqle  conducting  wedge.  Hutchins  and  Kouyoumjian  [16,17]  have  pre- 
sented a formula  for  the  diffracted  field  which  significantly  improves 
the  accuracy  of  the  solution  over  that  obtained  using  Pauli's  form.  This 
improved  solution  provides  superior  results  in  the  transition  regions  (near 
the  incident  and  reflected  shadow  boundaries). 

The  three-dimensional  wedge  diffraction  geometry  is  depicted  in  Fig- 
ure 3.  The  source,  located  at  s' (p1 ,4)' ,z' ) , generates  a radiated'?  field 
given  by  ^(s).  The  source  can  be  an  arbitrary  electric  or  magnetic  source 
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causing  cylindrical,  conical,  or  spherical  wave  incidence  on  the  wedge. 

The  diffracted  vector  field  at  s(p,<j),z)  can  be  written  in  terms  of  a dyadic 
diffraction  coefficient. 

Kouyoumjian  and  Pathak  [18]  have  derived  a more  rigorous  basis  for 
the  GTD  formulation  and  have  shown  that  the  diffracted  field  may  be  written 
compactly  if  it  is  defined  in  terms  of  a ray  fixed  coordinate  system.  This 
ray  fixed  coordinate  system  is  centered  at  the  point  of  diffraction  Q1, 

(or  points  of  diffraction  as  in  the  case  of  plane  wave  incidence).  Note 
that  Q'  is  a unique  point  or  set  of  points  for  a given  source  and  observation 
point.  The  incident  ray  diffracts  as  a cone  of  rays  such  that  the  cone 
half  angle  (rq)  equals  the  angle  (g^)  which  the  incident  ray  forms  with 
the  edge  as  shown  in  Figure  3. 

The  orthogonal  unit  vectors  associated  with  these  coordinates  (s', 

A A /A  A 

Bq.*'.s,b0i$)  are  related  by 

A A 

I = -S 


i = e;  x 

A A A 

s = Bo  X 4, 


where  I is  the  incident  direction  unit  vector  and  s is  the  unit  vector 
in  the  direction  of  diffraction. 


The  diffracted  field  is,  then,  given  by 


where 


+ 


In  addition,  N is  the  integer  which  comes  closest  to  satisfying  the  fol- 
lowing equations: 

?nnN+  - 8 = + tt 

?nuN  - B 3 -it 

with 

R = * - <|>  ’ 
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F(x)  is  defined  by 


F(x)  = 2j  x ejx  / e'j  T dT 

if*  I 

and  is  called  the  transition  function.  The  quantity  n is  defined  by  the 
wedge  angle  WA  where  WA  = (2-n)n.  The  (<|>-<t>')  terms  are  associated  with 
the  incident  field;  whereas,  the  C <j>+d> ’ ) terms  are  associated  with  the  re- 
flected field.  The  quantities  A(s)  and  L are  defined  later. 

In  matrix  notation  the  diffracted  field  is  expressed  by  [16,1 7] 


E?(sl 

-D$  0 

V(Q')  ' 

■eJ(s). 

• 0 -Dh  • 

.e|(Q') . 

(17) 


The  D$  coefficient  corresponds  to  the  E-field  component  parallel  to  the 
edge  with  the  (acoustically  soft)  boundary  condition 


Wedge  " °* 

The  Dh  coefficient  corresponds  to  the  E-field  component  perpendicular  to 
the  edge  with  the  (acoustically  hard)  boundary  conditions 


3E 
3 n 


= 0. 


Wedge 

The  quantity  A(s)  in  Equation  (17)  is  the  ray  divergent  factor  given 
in  general  by  fO] 


A<s>  -Jtfkr 
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For  the  straight  wedge  casep=p^,  where  p^  Is  the  radius  of  curvature 
of  the  Incident  wavefront  In  the  plane  containing  the  Incident  ray  and 
the  edge.  For  spherical  wave  incidence  A(s)  is  given  by 


A(s)  = 


The  quantity  L is  a distance  parameter  and  is  given  in  general  by 

Go] 

SVSV?  s,n?Bo 

L ' ce(0l*s)  (p,*s)  ('8 

where  and  p,,  are  the  principal  radii  of  curvature  of  the  wavefront. 
For  a straight  edge,  L is  given  by  [10] 


L _ s's  sin  8q 
s+s' 

for  spherical  wave  incidence. 


For  qrazing  incidence  (ip  * =0)  Ds=0,  and  must  be  halved  since  the 
incident  and  reflected  fields  merge  together  and  only  half  of  the  total 
field  on  the  surface  is  associated  with  the  incident  field  with  the  other 
half  being  the  reflected  field. 


0.  Curved  Wedge  Diffraction 


The  diffracted  field  due  to  a curved  wedge  is  analyzed  using  the  GTD 
techniques  developed  by  Kouyoumjian  and  Pathak  [IQ  . The  diffracted  field 
from  a curved  wedge,  as  seen  in  Figure  4 is  expressed  by 


id(d)-  i'(Q').  8 /infer  5‘J‘ 


1? 


OBSERVATION 

POINT 


(a) 


Figure  4.  Geometry  for  three-dimensional  curved 
wedge  diffraction. 


The  quantity  p is  the  caustic  distance  between  the  caustic  at  the  edge 
and  the  second  caustic  of  the  diffracted  ray.  This  can  be  found  using 


ne.(I-s) 

* — 

a sin  $ 
e o 


(20) 


A 

where  n0  is  the  unit  normal  to  the  edge  directed  away  from  the  center  of 
curvature  of  the  edge  at  the  diffraction  point  and  a > 0 is  the  radius 
of  curvature  of  the  edge  at  Q'.  The  quantity  pg  is  the  radius  of  curva- 
ture of  the  incident  wavefront  in  the  plane  of  the  incident  ray  and  tan- 
gent to  edge  at  Q1. 


The  diffraction  coefficients  for  the  curved  wedge  case  are  similar 
to  those  of  the  straight  wedge  in  Equation  (16).  The  only  major  modifi- 
cation necessary  for  the  curved  wedge  case  is  the  distance  parameter  L, 
which  appears  in  the  argument  of  the  transition  function.  The  diffraction 
coefficients  for  this  case  are  given  by 

2sin(J)  F kL1  a(d>-4> ' ) 
cos(^)  - cos(^-  ) 

t|c°t(j±|^jF[klrnaV4.')]  * 

+ j F[kLr0  aUV)]j  (21) 


DS(<M',80) 


4 


2nj2iik  sinB 


in  which  a(y)  = 2 coS^y  and  a+(y)  = 2cos^ 
distance  parameter  is  expressed  by 

Li  . s(Pe+s)  pJpJ  sin%o 

pe(p|+s)(p^+s) 


The  incident  field 


(22) 
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The  reflected  field  distance  parameters  from  the  surface  with  superscript 
o (reflection  boundary  at  tt-<|>')  and  from  the  surface  with  superscript  n 
(reflection  boundary  at  ( 2n-l ) are  given  by 

r s(Pg+s)  pjp£  s1n2BQ 

P<?(pl+S)  (p2+s)  ^ 

The  parameters  pj  and  p^  are  the  principal  radii  of  curvature  of  the  in- 
cident wavefront  at  the  point  of  diffraction.  For  plane,  cylindrical, 
and  conical  waves  p^  is  infinite  and  for  spherical  waves  Pg=s'.  The  pa- 
rameters pf  and  p£  are  the  principal  radii  of  curvature  of  the  reflected 
wavefront  at  O'  and  can  be  found  using  Equation  (7). 
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III.  CURRENT  DENSITY  INDUCED  UPON  A 
TWO  DIMENSIONAL  WAVEGUIDE 

A.  Introduction 


In  this  chapter  the  two  dimensional  waveguide  and  open  resonator  prob- 
lem is  investigated.  The  desired  goal  is  to  show  where  a simple  ray  pic- 
ture is  valid  when  compared  with  a more  general  modal  solution. 

8.  Waveguide  Mode  Analysis 

The  geometry  for  the  waveguide  mode  analysis  is  shown  in  Figure  5. 

In  this  section  the  radiated  field  from  the  parallel  plate  aperture  is 
obtained  by  superposition  of  the  diffracted  fields  from  the  edge  due  to 
each  propagating  mode.  Rudduck  and  Tsai  [19]  have  shown  that  this  method 
yields  answers  which  agree  with  a Wiener-Hopf  technique  and  experimental 
results.  The  waveguide  spacing  determines  the  number  of  propagating  modes 
which  exist  in  the  waveguide.  If  the  waveguide  spacing  is  "a"  then  the 
m-th  propagating  mode  has  a propagation  angle  (0  ) defined  by 


1 > s)n«„  - £ . (24) 

For  m 3 0,1,2 ^ , 

The  diffracted  field  for  one  mode  can  be  expressed  using  Equation 
f!5)  as 


i 

u = u 


z=t 


D 


77 


where  r * |r|  is  the  distance  to  the  far  field  observer  and  given  by 

Jk  ^ a 

r = | d | - I*.  • d|  and  i is  a vector  from  the  origin  to  the  point  of 
diffraction. 
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Figure  5.  Parallel  plate  geometry. 


The  parameter  D is  the  plane  wave  diffraction  coefficient  obtained  from 
Equation  (16).  The  incident  field  u1  can  be  written  as  the  sum  of  two  cros- 
sing plane  waves  in  the  waveguide.  For  the  transverse  magnetic  case 

, -jk(zcose  +xsinem)  - jk(zcos0  -xsine), 

H « £ [e  m m + e m m J 

, -jk(zcose  +xsinem)  -jk(zcos0  -xsine  ) 

Ex  ■ I + « " ■ ] 


E,  = z sinem  f-e 
z 2 o m *- 


-jk(zcose  + xsine)  -jk(zcosem  - xsine  ) 


At  the  edge  of  this  waveguide,  z=t  and  z=a, 
-jkt  cose_ 

V' 


Ex  ‘ V°%  e 


-jkt  cose. 


Ez  = 0. 


r 


1 


For  the  transverse  electric  case 


, -jk(zcos0m  + xsin0  ) -jk(zcos0m  - xsin0  ) 

Ey  = ± [-e  m m + e * " ] 

cos0m  -jk(zcos0m  + xsin0)  -jk(zcos0m  - x sin0  ) 

H = ”L  [e  m m - e m m ] 

x zo 

sin0„,  -jk(zcos0m  + xsin0_)  -jk(zcos0m  - xsin0l 
Hz  = —JL  [e  m m - e m m ] (27) 


and  at  z=t  and  x=a. 


Ey-0 


Hx  =0 


sine  -jkt  c0s8» 


H = 

z z 


m e 


(28) 


The  total  field  is  found  by  summing  in  the  far  field  the  individual 
diffraction  components  due  to  each  propagating  mode  in  the  waveguide.  Using 
reciprocity  this  will  yield  the  surface  current  or  charge  induced  upon 
the  waveguide  walls  due  to  an  incident  plane  wave  (see  Chapter  I). 

The  waveguide  mode  analysis  becomes  a numerical  accuracy  problem  when 
the  plate  spacing  is  large  in  terms  of  wavelengths.  As  the  plate  spacing 
is  increased  the  number  of  propagating  modes  increases  (see  Equation  (24)). 
For  large  plate  spacings  the  GTD  solution  described  in  the  next  section 
is  more  accurate  as  will  be  seen  later. 

C.  GTD  Parallel  Plate  Solution 

The  GTD  solution  for  the  parallel  plate  geometry  is  valid  when  this 
plate  spacing  is  large  in  terms  of  wavelengths.  The  total  electric  field 
is  given  by 
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where  m is  the  number  of  bounces  off  of  the  ground  plane  the  reflected  ray 
from  the  source  makes  before  reaching  the  far  field  and  n is  the  number 
of  bounces  which  occur  before  a reflected  ray  strikes  the  edge. 


The  incident  field  from  the  infinitesimal  source  located  in  the  wave- 
quide  is  shadowed  by  the  waveguide  geometry  for  all  far  field  observation 
directions  except  for  ep  = 0°  (see  Figure  F).  The  reflected  field  com- 
ponents are  found  using  image  theory  with  only  one  image  being  visible 
for  a given  range  of  ep.  Given  ep  the  desired  image  is  located  at 
-?a(n  - Tj-) x + O.z.  The  integer  n can  be  found  from 


2a(n-l) 

t 


< tan  fl  < 

P - 


2an 

t 


which  implies 
tan9 

n < — p-2  + 1 (29) 


where 


The  reflected  field  can  then  be  written  as 

0?m-l  -jk?amsinop  g-jkr 
K 6 


JV 


(30) 


where  R=1  for  H-field  parallel  to  the  edge  arid  R=-l  for  E-field  parallel 
to  the  edge. 


The  diffracted  field  is  composed  of  a sum  of  diffraction  terms  which 
are  generated  by  the  image  sources  as  in  Figure  6.  Given  by 
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Figure  6.  Imaqe  locations  for  the  parallel  plate  geometry. 
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where  £ is  a vector  from  the  origin  to  the  point  of  diffraction. 


D.  Results 

Figures  7 through  24  compare  the  GTD  current  solution  with  the  modal 
current  solution  for  various  waveguide  spacings,  a,  and  distances  from 
the  waveguide  aperture,  t,  as  defined  in  Figure  5.  The  parameters  a and 
t are  varied  independently  from  0.25X,  0.50X,  1.0X,  2. OX,  5. OX,  to  10. OX. 
Figures  7 through  15  are  for  the  H-f i el d parallel  to  the  edge  or  the  TM 
waveguide  mode,  while  Figures  16  through  24  are  for  the  H-field  perpendic- 
ular to  the  edge  or  the  TE  waveguide  mode.  Each  polar  plot  fixes  a,  t, 
and  the  polarization  and  varies  0 , the  angle  of  incidence  of  the  incoming 
unit  plane  wave,  is  varied  from  0°  to  90°. 


The  agreement  between  the  GTD  and  the  modal  solution  increases  as 
the  waveguide  spacing  increases.  As  the  spacing  is  increased  more  oropa- 
gating  modes  can  occur  in  the  waveguide  which  increases  the  accuracy  of 
the  GTD  rav  picture.  For  the  H-perpendicul ar  case  with  a=0.?5X  or  a=0.BX 
and  propagating  TE  modes  exist  in  the  waveguide  (see  Eguation  (?a))  so 
no  modal  solution  is  presented  in  these  plots.  The  effect  of  no  propa- 
gating modes  for  this  case  is  indicated  in  the  GTD  solution  by  the  field 
strength  being  very  low  for  these  cases. 


The  agreement  between  the  GTD  and  the  modal  solution  lessens  as  P^ 
approaches  90°.  This  effect  is  due  to  the  GTD  solution  approaching  a caustic 
since  the  number  of  ray  terms  is  approaching  infinity  as  can  be  seen  using 
Equation  (29)  with  0p  = 90°. 


Figures  25a  through  30a  show  the  z current  magnitude  in  dB  plotted 

i ~ 

for  the  H-parallel  case  normalized  to  the  incident  field  or  |J  |/|H  y|. 
Figures  25b  through  30b  show  the  y current  magnitude  in  dB  for  the  H-per- 


1 * i 

pendicular  case  normalized  to  the  incident  field  or  |J  |/|-H  cosP  x + H 

" - y p 

sinPpZ|.  On  each  plot  in  this  section  the  waveguide  spacing,  a,  is  fixed 
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and  the  distance  from  the  aperture,  t,  is  varied  continuously  from  0.3 
wavelengths  to  10.0  wavelenqths  for  four  different  values  of  plane  wave 
incidence  angles,  0p.  The  angles,  0 , occur  in  five  degree  increments 
starting  at  either  0°,  15°,  30°,  45°,  60°,  or  75°.  In  the  H-perpendicular 
case  for  a=0.?5  wavelengths  and  a=0.5  wavelengths  the  waveguide  is  in  the 
cutoff  region.  This  can  be  seen  in  the  plots  by  the  low  level  of  current 
indicated  in  Figures  43  and  44  and  by  the  exponential  decay  of  the  current 
away  from  the  aperture  as  seen  in  Figures  45  and  47.  As  stated  before 
the  current  value  presented  in  Figures  35  through  60  become  less  accurate 

as  0 approaches  90°. 
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Figure  A.  H„  in  the  waveguide  for  a=0.?5 
and  0.5  wavelength. 
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Figure  10.  H„  in  the  waveguide  for  a=I.O  wavelength. 
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Figure  12.  H„  i 

n the  waveguide  for  a=?.0  wavelength. 
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Figure  13.  H„  in  the  waveguide  for  a=5.0  wavelength. 
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Figure  14.  H„  in  the  waveguide  for  a=5.0  and  10.0 
wavelengths. 
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Figure  15.  H in  the  waveguide  for  a=10.0  wavelength. 
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Figure  16.  H in  the  waveguide  for  a=0.25  wavelength. 
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Figure  17.  ha  in  the  waveguide  for  a=0.?Fi  and  n.s 
wavelength. 
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Figure  1C».  H in  the  waveguide  for  a=1.0  wavelength. 
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Figure  22.  in  the  waveguide  for  a=5.0  wavelength. 
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Figure  23.  in  the  waveguide  for  a=5.0  and  10. n 
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Figure  24.  Hx  in  the  waveguide  for  a=l0.0  wavelength. 
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Figure  ?5b.  H±  current  density  for  parallel  plate  waveguide 
with  a=0.?5  wavelength. 
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IV.  CURRENT  AND  CHARGE  DENSITY  INDUCED  ON  A 
CYLINDER  ABOVE  A GROUND  PLANE 


A.  Introduction 


The  qeometrical  theory  of  diffraction  as  discussed  in  Chapter  II  is 
applied  to  the  cylindrical  geometry  illustrated  in  Figure  31.  The  restric- 
tions applied  to  this  model  are  that  the  source  must  be  located  at  least 
one  fourth  wavelength  away  from  the  cylinder  endcap  and  the  separation 
between  the  ground  plane  and  the  cylinder  must  be  large  enough  so  that 
the  waveguide  mode  nature  of  the  field  does  not  predominate  ( see  Chapter 
III).  For  patterns  in  the  x-z  plane  the  electric  field  is  given  by 

-»r  -‘r  r -*d 

E = E + EW  En 

n 

where  m is  the  number  of  bounces  off  of  the  ground  plane  as  the  energy 
propagates  from  the  source  to  the  far  field  receiver  and  n is  the  number 
of  bounces  which  occur  before  the  reflected  ray  strikes  the  edge. 

B.  Geometrical  Optics  Terms 


For  patterns  in  the  x-z  plane  of  Figure  31  the  incident  field  term 
is  shadowed  by  the  cylinder.  The  reflected  field  term  is  similar  to  the 
reflected  field  given  in  Chapter  II  in  that  Equations  f 3)  through  MA) 
are  applicable.  The  reflected  field  must  be  modified,  however,  to  take 
into  account  the  reduction  in  field  strength  due  to  reflections  off  of 
the  cylinder  which  defocuses  the  reflected  field  energy.  This  effect  can 
be  seen  using  Equation  (3) 

Er(s)  - Ei(Qr)  . U Ip.  e-J'ks  . 

> (pj+sHp^+s) 
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Figure  31 . Geometry  for  cylinder  over  ground  plane  case. 


In  general  the  effect  of  m bounces  off  of  the  cylinder  on  the  reflected 
field  is  given  by 


Er(s)  - e-Jk(s+2ms')  % 


/(D.o, 
Pi  +2 

p[(?W 

PlVlW  J( s+2msl)(s+pj,(m)) 


which  reduces  to 


Er(sl  = 


rm 

J*(2) 

pi 

°?. 

Y*  • • • 

p,nW  , 

(m-1 


- jks  ..  2 ? 

e'jkd,)l  (3?) 


p p 

in  the  far  field.  ■ The  quantities  and  ^ are  given  by 


rfll  2s'r. 
Pi  - 1 

1 ~ — nm 


rf£l  ^(Pi^'^^S1) 

P1  ~~  r for2<*<m. 

r1+2|p1(S'-1)+2s'  |sinep 

C.  Diffracted  Field 

The  diffracted  field  is  given  by  Equation  (19)  as 


In  the  above  equation  E 1 ( Q * ) is  the  reflected  field  incident  upon  the  edge 
of  the  cylinder.  The  parameter  p is  given  by  Equation  (20)  and  reduces 
to 

. _ P1  rl 

p _ - - . 

rj+p^sin^+sine  ) 

The  dyadic  diffraction  coefficient  is  given  by  Equation  (21)  with  L1 
= lrn  = Lr0  = p^.  For  the  principal  plane  patterns  taken  in  this  report 
it  is  sufficient  to  include  only  the  diffraction  from  the  point  on  the 
wedge  located  at  x=-r^,  z-t  in  regions  away  from  caustics. 

0.  Results 

The  procedure  described  in  this  chapter  assumes  that  the  diffraction 
off  the  endcap  is  dominated  by  the  principal  point  of  diffraction  located 
at  x=r^ , y=0.,  z=t.  This  assumption  is  usually  valid  if  the  other  points 
of  diffraction  are  separated  from  the  principal  diffraction  point.  If 
these  diffraction  points  are  sufficiently  close  together  a caustic  field 
can  occur  yielding  an  erroneous  result.  These  caustics  become  more  of 
a problem  as  the  cylinder  radius  is  decreased.  The  results  for  this  sec- 
tion are  included  with  the  results  given  in  Chapter  V.  Chapter  V describes 
a technique  which  can  be  used  to  partially  correct  these  caustics  caused 
by  closely  spaced  diffraction  points  around  the  endcap  of  the  cylinder. 


V.  CAUSTIC  CORRECTION 


A.  Introduction 

Due  to  the  large  number  of  rays  necessary  to  obtain  a GTD  solution 
to  the  cylinder  mode,  there  is  a good  possibility  of  introducing  caustics 
into  the  pattern.  These  caustics  occur  when  three  diffraction  points  on 
the  endcap  become  close  to  each  other.  These  caustics  can  be  partially 
corrected  using  the  method  of  Albertson,  Balling,  and  Jenson  [20] . 

B.  Theoretical  Section 
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Figure  32.  Geometry  for  the  caustic  correction  procedure. 
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Substituting  this  expression  into  the  phase  integral 
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The  function  D ^|t  e ^is  the  parabolic  cylinder  fund 
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proximated  by  [23] 

D j (t  eJ  4 ) = 1 . 215280214y1( t)  - 0.581368317y?(t) 
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s the  parabolic  cylinder  function  and  is  ap- 
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The  maximum  possible  error  due  to  including  up  to  the  eighth  order  term 
in  the  parabolic  cylinder  function  is  less  than  1%  for  |t|<2. 

The  above  result  may  be  expressed  as  a caustic  correction  factor, 

c,  used  to  multiply  the  GTO  solution  in  a case  caustic  region.  In  a region 

where  there  are  four  diffraction  points  and  h4(<K)  > 0 

, 3n  z 
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with  a = hM(4>2)  V 3k/h^4'(<t»^).  In  the  region  where  there  are  two  diffraction 
points  and  h^  (<f>£)  > 0 
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For  i/4^^)  < 0 the  complex  conjugate  of  the  correction  factors  should 

he  used.  A distinction  between  central  and  noncentral  rays  must  also  be 

made  in  the  case  where  four  rays  are  present  (see  Figure  33).  If  the  ray 

is  noncentral  then  h"(*i)  is  replaced  by  -2  h"(<j)t)  in  evaluating  o.  Non- 

' (41  c 

central  rays  can  be  recognized  if  h"(<j>^)  * hv  ' ($,!,)  > 0.  For  the  central 
ray  if  h"(<t>Jp  * f/4^  (<(>£)  < 0 four  rays  are  present  and  if  h"(4>£)  * (/4^ 

(<fp)  > 0 two  rays  are  present.  Note  that  the  parameter  a can  he  used  to 
check  if  a caustic  field  is  present.  If  o<2  then  a caustic  effect  is  present 
and  if  o> 2 then  the  caustic  effect  is  negligible  and  c may  be  approximated 
by  one. 

For  the  geometry  illustrated  in  Figure  22 
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T =Jxx+^s+zs  +ri'^r1xscos4) 1 - 2r1yssin4>1  . 

C.  Results 

Figure  34  and  35  show  the  effect  of  the  caustic  correction  on  a five 
wavelength  radius  cylinder  located  one  wavelength  above  the  ground  plane. 
Since  only  the  principal  point  of  diffraction  is  used  in  the  caustic  cor- 
rection, the  current  level  in  the  caustic  region  is  slightly  below  the 
level  expected  in  the  region  where  three  diffraction  points  are  present. 
The  effects  of  the  other  two  diffraction  points  cannot  he  included  since 
they  would  he  diffracted  rays  which  are  launched  as  creepinq  waves  along 
the  cylinder.  A solution  for  this  type  of  ray  is  not  currently  available. 

Figures  35  through  71  follow  the  same  format  as  Figures  ?5  through 
30  in  presentinq  the  induced  current  on  the  cylinder,  except  that  the  cy- 
linder radius  is  varied  from  0.?5X,  0.5X,  1.0X,  ?.0X,  5.0X,  to  10. OX. 
Figures  36a  through  71a  show  the  axial  current  density  in  dB  normalized 
to  the  incident  plane  wave  or  |Jz|/|H1y|.  In  Figures  36b  through  71b  the 

surface  charge  density  is  plotted  normalized  to  the  incident  plane  wave 
i A i a 

or  |p|/|eE  cos0pX-eE  s i n© pz | . In  Figures  36c  through  71c  the  circumfer- 
ential current  density  is  presented  normalized  to  the  incident  plane  wave 
i ^ i A 

or  | J | / | - H cos0pX+H  sin0pZ|.  The  caustic  corrected  solution  is  used  for 
all  these  computations  in  that  it  provides  the  most  accurate  results. 

The  caustic  problems  become  more  pronounced  as  the  cylinder  radius  becomes 
smaller.  When  the  radius  is  ten  wavelengths  the  solution  is  very  similar 
to  the  parallel  plate  waveguide  as  expected.  The  current  solution  becomes 
less  accurate  as  the  plane  wave  incident  angle,  0p,  approaches  90°  similar 
to  the  parallel  plate  waveuide  solution.  This  effect  is  reduced  somewhat 
due  to  the  reduced  enerqv  of  the  ravs  after  bouncing  several  times  off 
of  the  cylinder  surface. 
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Figure  40c.  Circumferential  current  density  with  a=0.?5  wavelength 

and  Tj/ a=?0.0. 
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Figure  42 b.  Surface  charge  density  with  a=0.5  wavelength 

and  r./aO.S. 
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Figure  44a.  Axial  current  density  with  a=0.5  wavelength 

and  r1/a=?.0. 
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Figure  47h.  Surface  charge  density  with  a=0.*  wavelength 

and  r1/a=?0.0. 
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Figure  50c.  Circumferential  current  density  with  a=1.0  wavelength 

and  r. /a=l .0. 
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Figure  55a.  Axial  current  density  with  a=2.0  wavelength 

and  rj/a=0.25. 
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Figure  56c.  Circumferential  current  density  with  a=2.0  wavelength 

and  r./a=0.5. 
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Figure  57c.  Circumferential  current  density  with  a=?.Q  wavelength 

and  r1/a=1.0. 
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Figure  71c.  Circumferential  current  density  with  a=10.0  wavelenqth 

and  rj/a=1.0. 


VI.  CONCLUSIONS 


The  object  of  this  research  has  been  to  analyze  the  surface  current 
and  charge  density  induced  on  a perfectly  conducting  cylinder  situated 
above  a perfectly  conducting  ground  plane  by  an  incident  plane  wave.  This 
model  can  be  used  to  simulate  a parked  aircraft  on  a runway  or  other  similar 
structures.  The  basic  approach  used  in  this  study  has  been  to  apply  the 
Geometrical  Theory  of  Diffraction  to  analyze  this  problem.  Thus,  the  solu- 
tion presented  here  for  surface  current  and  charge  density  induced  on  a 
finite  cylinder  above  a ground  plane  is  applicable  for  radii  as  small  as 
one  quarter  wavelength. 

In  Chapter  III  the  parallel  plate  aperture  problem  was  investigated 
in  order  to  indicate  where  the  GTD  ray  picture  was  valid  and  to  verify 
the  accuracy  of  the  solution.  The  GTD  results  compare  very  favorably  with 
a modal  type  solution.  Current  densities  for  a variety  of  aperture  spacings, 
and  current  locations  along  the  waveguide  walls  were  presented. 

In  Chapter  IV  the  results  for  the  cylindrical  geometry  were  presented. 

In  this  section  the  charge  and  current  density  were  plotted  as  a function 
of  the  cylinder's  distance  above  the  ground  plane,  the  cylinder  radius, 
the  incoming  plane  wave's  angle  of  incidence  with  respect  to  the  ground 
plane,  and  the  distance  from  the  observation  point  to  the  cylinder  endcap 
along  the  bottom  of  the  cylinder.  In  general  the  results  have  shown  that 
caustics  can  appear  in  the  solution  especially  if  the  cylinder  radius  is 
small  in  terms  of  wavelengths.  A partial  correction  for  these  caustics 
was  discussed  in  Chapter  V and  is  included  in  the  results.  The  caustic 
correction  procedure  is  not  complete  since  using  the  current  theory  only 
the  principal  diffraction  point  can  be  used  in  the  procedure. 


Incident  plane  waves  in  the  principal  plane  were  used  in  this  study. 
However,  the  theory  can  be  extended  to  include  incident  plane  waves  out 
of  the  principal  plane.  There  is  also  a need  for  experimental  data  in 
this  area  to  verify  results  obtained. 
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